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I. INTRODUCTION 

Quantum theories whose evolution was driven by a 
non-Hermitian Hamiltonian have a long history emerging 
in the early times of quantum physics (for the earliest ac- 
counts see [l| where the indefinite metric quantum theo- 
ries are discussed in detail; modern approaches are based 
on treatment of 0). They became a prominent topic 
of discussions in last decade after publication of the key 
article Q mostly under the name "PT-symmetric quan- 
tum theories. A vital point of modern interpretations of 
such theories is always the redefinition of the scalar prod- 
uct which, assuming real spectrum of the Hamiltonian, 
enables to re-establish the consistent probabilistic inter- 
pretation of the theory. The standard machinery goes as 
follows: At first, one takes the theory's non-Hermitian 
Hamiltonian H and check whether it has real spectrum. 
If so, one tries to find an operator O, called a metric, 
which satisfies 

H^Q = QH. (1) 

The metric has to be positive, bounded, invertible with 
bounded inverse, Hermitian operator, and must map the 
range of H into the domain of definition of . (For 
Hamiltonians with complex spectrum such metric does 
not exist.) Having found one may define new scalar 
product 

(*|$) e = (*|9$> (2) 

which should be used for calculating transition ampli- 
tudes. Equivalently, if one is able to decompose the met- 
ric into 

9 = fttfi, (3) 

one can define a new, Hermitian Hamiltonian by a simi- 
larity transformation as 

h = QHQ- 1 . (4) 



'Electronic address: bila@ujf.cas.cz 



If (HJ) holds for some operator A in place of H, we say 
that A is quasi-Hermitian with respect to 0. All ob- 
servables of the theory have to be quasi-Hermitian with 
respect to the same metric. For any quasi-Hermitian ob- 
servable we can of course define its Hermitian counterpart 
a = QAQ^ 1 . Although the Hermitian representation us- 
ing ftasa Hamiltonian is assumed to exist for any model, 
for practical purposes it is often preferable to make cal- 
culations in the original basis using the metric, since the 
h (and eventually other observables) tend to be rather 
complicated and inconvenient to deal with. 

It is not difficult to see that the probabilities are con- 
served under described setting, for we have 

— (#|0$) =i(HV\G$) -i(*|9ff$) (5) 

which vanishes due to quasi-Hermiticity condition. The 
quasi- Hermiticity condition can be rewritten in terms of 
the evolution operator 

U(t',t) = e~ iH(t '- t \ (6) 

Instead of standard U'U — I one has now 

U^OU = Q. (7) 

It may be noted that the equation (fT]) does not de- 
fine the metric uniquely. To specify a unique metric, one 
has to employ the quasi-Hermiticity condition for other 
observables; if the set of observables is irreducible there 
can be only one metric which fulfills the given require- 
ments. The situation can become more peculiar when 
the Hamiltonian is time-dependent. In such a case there 
is a continuous set (depending on t) of equations 

H(tfe = eH(t) (8) 

that have to be satisfied. The set {if (t)} can well be irre- 
ducible and specify the metric uniquely, but it is entirely 
possible for the metric to be over-specified and thus non- 
existent. This leads to the concept of quasi-stationarity: 
the Hamiltonian is said to be quasi-stationary if there 
exists a single metric with respect to which it is quasi- 
Hermitian at any t. Only quasi-stationary Hamiltonians 
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are accessible for the standard physical interpretations. 
Note that one cannot simply assume a time-dependent 
metric since it would lead to an additional term in ([5]) 
proportional to dO/di. We will discuss how to overcome 
these problems in section Hill 



II. SCATTERING IN NON-HERMITIAN 
THEORIES 

The number of 'PT-symmetric models investigated in 
recent literature is vast. Most of the research has been 
concerned mainly with spectral properties of the systems 
and especially the presence of bound states. Whilst the 
treatment of bound states has been standardised and 
became a routine task, the discussion of scattering has 
faced several conceptual difficulties. The problems are 
closely related to the metric redefinition, of course. To 
see the core of the complications, let us consider two spe- 
cial cases: 

At first we look more closely to quantum-mechanical 
scattering with short-range potential. The standard ap- 
proach to the problem uses x-representation, solving the 
Schrodinger equation with plain- wave-like boundary con- 
ditions for specification of the incoming wave at large 
\x\ in a chosen direction; the scattering amplitudes are 
taken as the asymptotics of the solution again at large 
\x\. This is reasonable since in asymptotic regions (out of 
range of the interaction) the solutions from the continu- 
ous spectrum solve the free-particle Schrodinger equation 
and thus can be interpreted as free particles coming into 
(or out from) the interaction region. 

Such approach can hardly be justified with non- 
Hermitian interaction. The reason is that x in the wave 
function ij)(x) has no longer the interpretation of an 
eigenvalue of the position operator. This is mainly due 
to the scalar product redefinition: the operator X acting 
as [Xi/j](x) = xtp(x) is generically not 0-quasi-Hermitian 
and therefore not acceptable as an observable. In order 
to find the real asymptotic states one has to construct 
physical position operator first. This can prove difficult 
and even impossible in some situations, namely within 
relativistic quantum mechanics [l2l ]. 

When the use of ^-representation is inconvenient 
(which includes most of the field theory, relativistic and 
many-particle physics), it became customary to use dif- 
ferent approach, which is based on decomposition of the 
Hamiltonian into free part Ho and interaction part H\. 
The scattering amplitudes are the matrix elements of the 
S-matrix between the free Hamiltonian eigenstates; S- 
matrix is given by 

S= lim U (0,t + )U(t+,t-)U (t-,0) (9) 

t± — >±oc 

where the indices differentiate between free and full evo- 
lution operators in the same way as for the Hamilto- 
nian. The sense of the convergence in © has to be 
clarified. For physical reasons we can expect the limit 



to exist for particular matrix elements. In many-particle 
systems these are the matrix elements between particle 
states whose momentum-representation wave functions 
have disjoint supports at t = 0; this is because for large 
\t\ the free evolution operator removes them far apart 
and hence out from the range of interaction; after that 
they evolve exactly the same way under Uq as under U 
and the contributions of large \t\ on the right-hand side 
of ([9]) cancel. 

With a non-Hermitian Hi the situation is more compli- 
cated due to the necessity of scalar product redefinition. 
The main problem is usually the fact that H and Hq, 
although both quasi-Hcrmitian, are not so with respect 
to the same metric. This represents a serious obstacle 
because ©, however formally treated, would lead to an 
S-matrix unitary with respect neither to Oo nor to O (0o 
is the free metric, i.e. calculated for the free Hamiltonian, 
and often it is simply unit operator). 

A useful trick which helps to get better meaning to 
© is the introduction of adiabatic switching of the in- 
teraction. This consists of introduction of an artificial 
time-dependent damping factor to the interaction Hamil- 
tonian. The exact form of this factor is not important as 
long as its derivative is small and the factor is equal to 1 
at t = and tends to zero at infinities. For concreteness 
one can take 

H s (t)=H +e- s ^H l (10) 

with e small positive. Hence for large |t| the full and free 
Hamiltonians are essentialy the same and the limit © 
may exist on larger domain [131 ]. The change presum- 
ably does not affect the scattering amplitudes between 
particles that are outside the range of interaction when 
e _£ ''! differs significantly from one. The original scatter- 
ing problem can be recovered by putting e — > ■ 

Analogous approach to the problem of scattering for 
non-Hermitian H would face all the problems of time- 
dependent Hamiltonians described in the introduction. 
Since in general there is little chance to find common 
metric for H and H n (and thus for all H E (t) producing the 
quasi-stationary scenario) we have to develop physical 
interpretation for systems with time-dependent metrics. 

III. TIME-DEPENDENT METRICS 

It turns out that in order to evade the demand of quasi- 
stationarity one cannot maintain H both as an observ- 
able and as the time evolution generator. In the latter 
case one finds the metric by solving ([TJ) for each t sep- 
arately and then postulates the Schrodinger equation in 
the Hermitian representation 

= mm (ii) 

(we keep the distinction between small and capital let- 
ters, which distinguishes the original and the Hermitian 
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representation, for the wave functions too, i.e. if — f^)- 
One must now sacrifice the validity of 



(12) 



since for time-dependent 6 the equations (fTTj) and (fT2|) 
are not equivalent. Instead of H, the time evolution in 
the non-Hcrmitian representation is now driven by 



Hgen = H — if} ^Cl 



(13) 



where dot denotes the time derivative. This resolution 
was suggested in 

Alternatively, one can maintain (|12[) and postulate the 
probability conservation 



dt 



;*|e$) = o 



(14) 



for all \P and <1>. Combining (fT2")) and (fT4"|) one gets 
© = i(QH - H^O). 



(15) 



Hence, the algebraic equation (TT|) is replaced by differen- 
tial equation (fTSf . The condition (fT5|) is needed only for 
the Hamiltonian, the observables must still satisfy ordi- 
nary quasi-Hermiticity condition. Consequently, if there 
is not a static (time independent) solution to ([15)) , H is 
no more an observable. However, we can construct an 
"observable Hamiltonian" 



H, 



obs 



H 



(16) 



If necessary, an equivalent Hermitian representation can 
be obtained. It has no sense to transform H, rather one 
shall apply the similarity transformation to the evolution 
operator U which yields 



u(t,t') = Q(t)U(t,t')Q(t'y 

Since we have 

h = i du(t,t>) x 
dt y ' 

we can easily realise that 

h = r^-f/obs^ . 



(17) 



(18) 



(19) 



The transition probability amplitudes are in this 
framework computed as ©(i)-induced scalar product of 
concerned Schrodinger-picture state vectors taken in the 
same time instant as the metric. For instance the prob- 
ability of finding state <I> at time t in a system prepared 
in state ^ at time t' is 

p = |($|e(t)[/(M')l*)l 2 - \{<s>\u(t',tfe{t')\y)\ 2 . 

(20) 

Let us note that the non-observable nature of H is in 
fact not a grave problem. When the Hamiltonian is time- 
dependent, the energy is no more conserved and cannot 



be directly measured even one we deals with completely 
conventional Hermitian system. The energy such a sys- 
tem has well-defined approximative meaning only if the 
time dependence is weak. In that case one can as well 
use -Hobs provided is small. That such slowly changing 
© exists is not a priori clear and we will return to this 
point in later discussion. 

The former approach with the observable H is suit- 
able for situations where one intends to use the non- 
Hermitian Hamiltonian as an auxiliary representation of 
a time-dependent Hermitian h which is given as the pri- 
mary operator. If the theory is formulated with H as 
a generator of time evolution (as it clearly is in case of 
the adiabatic interaction switching) one must accept the 
latter resolution, however it is more complicated. Com- 
pared treatments are not so different as it could appear 
on the first sight: The similarity of formulae (|13p and 
(Unj), as well as the validity of (fH?)) suggest that by iden- 
tifying H (of the former approach) with -ff bs and H gcn 
with H (of the latter approach) an isomorphism between 
these two methods can be established; it only matters 
where one starts and what has to be calculated. For rea- 
sons given above we will exclusively treat H as a time 
evolution generator in the following. 



Calculation of the metric 

Obviously the calculation of the metric based on (fTT 
is complicated and has to be done numerically or approx- 
imatively in majority of conceivable models. The appar- 
ent similarity between (|15[) and the Schrodinger equation 
for the evolution operator or the Heisenberg equations for 
observables suggests we can use similar techniques, e.g. 
iteration of the equation rewritten as 

0(t) = 0(O) + / Q(t')H(t')- H(t')^e{t')dt'. (21) 
Jo 

Contrary to the Schrodinger equation the process is fur- 
ther complicated by the ambiguity in setting O(0), which 
can be resolved only with more detailed knowledge of the 
specific model. However, if O(0) is chosen Hermitian, so 
is 0(i) due to the right-hand side of (fT5)) . Of course 
there is a possibility not to calculate O directly and find 
the evolution operator instead. The metric is then, as 
suggested by ([2U1) . 



e(t) = u{t',tye(t')u(t',t) 



(22) 



Although it is uneasy to give details about the metric 
in general case, hopefully we can say more about systems 
with slowly changing H{t) which are the main topic of 
interest in this paper. Even more specific case is that 
of time-independent Hamiltonian: one may ask about 
the solutions of (fT5|) when H does not depend on time. 
If H is quasi-Hermitian, which we will always suppose 
even if we can consistently solve (TT51) for broader class 
of Hamiltonians, then there is always, by definition of 
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quasi- Hermiticity, a static solution which we will denote 
s . But in addition there are other, non-static solutions. 
In finite-dimensional Hilbert space, there is always a basis 
in which their matrix elements are periodic. Indeed, if 
we write the metric as 



(23) 



where \f r ™ is the n-th eigenvector of W and $ mn = i?£ m , 
([T5"]) gives 



(24) 



which leads to constant diagonal elements and periodic 
off-diagonal elements. We have used the biorthonormal- 
ity relations between *S? n and the eigenvectors of H which 
we denote : 



(* m |*n> 



(25) 



The whole metric is periodic if there is a common multi- 
ple of all energy differences. Note that if H loses its real 
spectrum z? m „ will exponentially rise or decay (diagonal 
elements included) and the periodicity will be lost. Un- 
der specific conditions the argument can be extended to 
infinite-dimensional diagonalisable Hamiltonians. 

For constant Hamiltonian one can formally solve the 
equation ([2T]) . In special case with 8(0) — I one gets 



e(t) = I + i(H-H^)t- (H 2 -2H t H + H 



.2 J 2 



(26) 



the colon in the last term standing for normal ordering 
of the exponential's Taylor series, which means all 
standing left of any H in each monomial. With differ- 
ent starting metric one has to insert 0(0) in each term 
between the rightmost W and the leftmost H. Tak- 
ing truncated series (|26p makes a good approximation if, 
roughly speaking, (H — H^)t is small. 



Time-dependent metrics in scattering problems 

In order to make sense of scattering in theory with adi- 
abatically switched non-Hermitian interaction we need to 
modify the definition of S-matrix in @ to reflect the pres- 
ence of the metric. It will be convenient to separate the 
S-matrix into product of the M0ller operators 



lim U (0,t)U(t,0), 



= lim U{0,t)U o (t,0) 



(27a) 
(27b) 



and introduce the customary notation of in- and out- 
statcs as 



vj/ ln 

^out 



(28a) 
(28b) 



Then, the transition amplitudes are given as ^I^P) 
and may be easily generalised for use in non-Hermitian 
theories to 



s fi - w'!©!*! 11 ), 



(29) 



subscripts f and j denoting final and initial state as usu- 
ally. 

A convenient interpretation of in-states (out-states 
analogously) is that, under limit t— — » — oo they are 
states evolved first from to £_ under free Hamiltonian 
(which is motivated by need to get rid of the oscillating 
phase in [/(£_, 0) in the limit) and then back to t = 
under full Hamiltonian. Applying the logic of adiabatic 
switching we shall set the metric at t — — oo to conform 
the free Hamiltonian and let it evolve using (fTS"]) . Hence 
in (|29[) is not an arbitrary metric; on the contrary it 
is uniquely specified as the metric which adiabatically 
evolves from the free metric ©o (which we will often as- 
sume to be equal to I): 







lim lim U E (t,0)^e o U E (t,0), 

e—*0 t—> — oc 



(30) 



e having the same meaning as in (|10p . Note that U e is 
not unitary and so the right-hand side does not simply 
evaluate to / even in case of 0o = /. 

To make the approach consistent, we have to show that 
the limit in (|30[) exists and is unique. In particular this 
implies the following requests: 

1. When transition between two time- independent 
Hamiltonians happens adiabatically, a static met- 
ric must evolve into another static metric. E.g. 
changing Hq to H\ during some interval (0, T) with 

Hr{t) — ^T^- Ho + ifHi in the meantime, the met- 
ric which started as static for Hq must evolve into 
static metric for Hi when T — > oo (the last condi- 
tion replaces e — > here). 

2. The resulting metric should be independent on the 
precise manner of the switching, i.e. all different 
damping functions should lead to identical results 
in the adiabatic limit. 

We have seen that for time-independent quasi- 
Hermitian diagonalisable Hamiltonian the solution to 
(fTS")) is bounded for all t and it oscillates around the static 
solution s . In the (non-orthogonal) eigenbasis of 
the static solutions are diagonal and the non-diagonal 
elements are sinusoidal. The static solutions are there- 
fore not unstable, i.e. if 0(0) is close to any of the static 
solutions it remains close for all t. This is, of course, ben- 
eficial: otherwise the smallest deviation from the static 
metric, which necessarily occurs for finitely slow transi- 
tion, would lead to utter collapse of the staticity of the 
metric at later time, and to cure that in the adiabatic 
limit (i.e. infinitely slow transition) could be impossible. 

Our demands resemble somewhat the adiabatic theo- 
rem of ordinary quantum mechanics. Recall that the adi- 
abatic theorem says that the eigenvector of Hq evolves 
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into eigenvector of H\ if the transition is adiabatic. 
Standard proofs use the criterion of adiabaticity which, 
roughly speaking, says that dH/dt <C 6 during the whole 
transition, where S is a spectral gap separating the re- 
spective eigenvector from the rest of the spectrum. We 
can translate the theorem to suit our purposes, replac- 
ing "eigenvector" with "static metric" , albeit there is an 
important difference - the absence of gaps in the set of 
static metrics. 

Although we cannot generalise the proof of adiabatic 
theorem directly, we can indeed use it to establish the 
adiabaticity of the metric, even only for Hamiltonians 
with discrete non-degenerate spectrum. Moreover we will 
assume that for both Hq and H exists some 

6 = imn\E n -E m \ (31) 

and similar minimum exists for the transitional Hamil- 
tonian H e for all t. Then the gap condition holds for 
sufficiently slow transition and in the adiabatic limit we 
have 

limE/ E (-oc,0)* n = $„, (32) 

for all which are eigenvectors of Ho while $ n are 
eigenvectors of H . Thus due to ([23)1 we can write 

C/ E (-oo,0)=^|$ n )(*"|+O( £ ). (33) 

n 

A static metric for Hq can be written as 

e = X>»l $ ">( $n l ( 34 ) 

n 

(when Hq is Hcrmitian possibly <&" = <&„ and with 0„ — 
1 the unit metric is recovered) . Combining (|3~4")) with ([BIT]) 
using (|2"5j) yields 

e = J2$n\y n )(V n \+0(e), (35) 

n 

i.e. for e — > a metric which is static because diagonal 
in the eigenbasis of Note that (|35|) does not specify 
the metric uniquely since the biorthonormal eigenbasis 
can have several distinct normalisations and the normal- 
isation is not specified in (I32)) . 

IV. TOY MODELS 

Two-level system 

To illustrate the behaviour of time-dependent metrics 
we present two different toy models. The first one is a 
two-level system whose Hamiltonian will be parametrised 
by Pauli matrices: 

H = ho + hidi, (36a) 
9 = 0o + 0^, (36b) 



with % £ {1..3}. Equation (fl5| then reads 

0o + dm = 2i[(0 o + 0*0*) Im ho + (0i + ^qct<) Im h t ] 
+2a i e: ijk '& j Rehk, (37) 
e being the Levi-Civita symbol. After putting 

2Reh (i = Vf, (38a) 
21mh li = (38b) 

with v,w £ M and £ {0..3} we get the following equa- 
tions 

-0o = ^Wft (39a) 

-fy = ■d w i + -&iW + £ ljk tfjVk- (39b) 

Absolutely simplest is the case of Hcrmitian, time- 
independent H, i.e. — 0. The equations then reduce 
to 

= (40a) 

d = vxS. (40b) 

Here the 0th component decouples from the rest and can 
be adjusted to an arbitrary constant. The vector pre- 
cesses with constant velocity around the direction speci- 
fied by the vector v. If 0(0) and v are collinear, then the 
solution is constant - in such case H and commute. If 
not, the frequency of precession is proportional to v, in 
other words, to the difference between the energies, as it 
was expected from (IM)) . and it does not depend on the 
initial condition. 

The case u> M ^ (i.e. non-Hermitian, but still constant 
H) is analogous. Either it can be transformed by a sim- 
ilarity transformation to the previous case, or managed 
directly: First we find the static solution s , if possible. 
The static solution exists if H is pseudo-Hermitian, i.e. if 

wo = 0, (41a) 
v-w = 0. (41b) 

Then, 

s = --£[vx w] +av, (42) 
v 

where 0g and a can be chosen arbitrarily. General solu- 
tion is easily obtained by setting up the Ansatz 

= av + 0w + j[v x w]. (43) 

It follows that a — const, and, up to a time shift, 

(3 oc sin \/ v 2 — w 2 t. (44) 

The remaining unknowns we get from 7 = and 0o = 
—w 2 /3. This is again what we have expected: if v 2 > w 2 
the energies of H are complex and consequently the pe- 
riodicity of the metric is replaced by exponential growth. 

When the time dependence is enabled we again see 
the expected behaviour. When the Hamiltonian changes 
rapidly, the metric is disrupted from its static state and 
starts to oscillate, whereas under slow transition the os- 
cillations are suppressed. An illustration is given if Fig. [T] 




FIG. 1: Illustration of the evolution of metric of the discussed 2-dimensional system. For concreteness w = (0,0,3) is kept 
fixed as well as vq and V3, the remaining components of v are vi(t) = v%{T — t) = 2t/T between t = and t = T (the latter 
marked by the dashed vertical line in the graphs) and remain constant before and after. The graphs show the evolution of $0 
and three components of $, starting from the static solution i? s . If the change of H is fast, as shown in the left graph where 
T = 1, the metric starts to oscillate. On the other hand, an adiabatic change in the left plot (T = 100 is much greater than 
the period of the metric of time-independent system) leaves the metric almost in a static state. 



Cubic oscillator 

The second model we will ivestigate here is the har- 
monic oscillator with cubic perturbation on the real line, 



H = p 



159 



(45) 



with real g; p and q are standard position and momentum 
operators. This Hamiltonian belongs between the most 
researched non-Hermitian Hamiltonians (except 0, Q see 
also older |6|). It is unfortunately not exactly solvable 
(albeit there is a proof of the spectrum being real Q) 
but the metric can be calculated perturbatively. It has 
been done in || and we will in principle follow the way 
outlined there which treats the metric as a function of 
the operators p and q, although for sake of simplicity 
we will not pursue the most efficient evaluation method 
that parametrises the metric as an exponential of other 
operator. Instead we will write the metric directly as 



e = Q(p,q), 



(46) 



To resolve the problems of operator ordering the parame- 
ters p and q in (|46p and (145|) will be considered c-numbers 
during the calculations, using the Moyal product in place 
of operator multiplication. (We use the following defini- 
tion 



lF*G]( P ,q) = F(p,q)exp- 



d q dp 



dpd q 



G(p,q). (47) 



for the Moyal product; arrows show the direction in which 
the derivative acts. Its particular advantage is that if an 
operator A is represented by function A(p, q) then A^ is 
represented by conjugated A(p, q)*.) 

Prior to the perturbative evaluation of we can look 
at g = case which is the ordinary harmonic oscillator. 



Here the time-dependent metrics are obtainable exactly. 
In chosen representation (|15p is properly rewritten as 



e = i(e *h - ir *e) = 2{ q 



dp dq' 



(48) 



It seems favourable to move to polar coordinates in p-q 
plane, so that p — psinip and q = p cos ip. Then the 
p-derivatives vanish and (|35]) becomes 



~dt 



dp 

The general solution to (j49|) can be sought in form 
9(p,tp,t)=A(p)S(t + tp/2). 



(49) 



(50) 



Such function satisfies (|4"9"| when 2 is periodic with pe- 
riod 7r. Moreover, both A and 5 must be analytic if we are 
to be able to convert them into series of operator products 
and A shall be bounded in order to get bounded metric 
(we can apply stronger criteria for A to ensure that the 
metric is well behaved, equation (j49| puts no restrictions 
to the p-dependence of 0). The static solution here is 
any function s which depends only on p; such function 
corresponds to an operator which commutes with H for 
H(p,q) = p 2 . The periodicity of the harmonic-oscillator 
metric is clearly result of the equidistance of its spectrum. 

For non-zero g one has to rely on perturbative calcu- 
lation. Due to discussion in Q a static solution can be 
obtained by putting 



= l + g0 1 {p,q)+g 2 6 2 (p,q) + 



(51) 



with 0j being the most general polynomial of order i + 2. 
Thus, in the first order in g we get 

61 = a + {]p+~fq+5p 2 + epq + (q 2 + r/p 3 + tfp 2 q + Lpq 2 + nq 3 . 

(52) 
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Substituting it into (|15|) after a rather tedious calculation 
we get the first order contribution to the static solution, 
which is 



CONCLUDING REMARKS 



ef = c + d(p 2 + q 2 )+pq 2 + -p 3 , (53) 

c and d are arbitrary real numbers. Their presence re- 
flects the ambiguity of the metric. In fact, the first two 
terms correspond to an operator which commutes with 
Hq = p 2 + q 2 and hence they are contributions to the free 
metric 0o- Since we wish to start with Go = / we expect 
that the adiabatically evolved metric will have c = d = 0. 

Like in case of our 2-dimensional model, we have cal- 
culated the time-dependent metric for linear interaction 
switching, i.e. 



Hi = 
Hi = 



for t < 
itgq 3 



T 



for t e (0, T) . 



Hi = igq 3 for t > T. 
For initial condition 6(p, q) = 1 we arrive at 



(54a) 
(54b) 
(54c) 



e 



T 



— (24t - 27 sin t + sin 3t)p 3 
36 



■-(2 + cost) sm 4 -pV 



• sin t 



1 

12 



sin it)pq 



(55) 



--(15 + 2 cost + cos2i) sin 2 ^<? 3 



0(.g 2 



for £ G (0, T). Looking at ([55]) in < = T one can realise 
that in the adiabatic limit T oo only two terms sur- 
vive; they are equal to 1 + g<df + 0(g 2 ) with zero c and 
d as expected. 



We tried to argue that the scattering in 'PT-symmetric 
quantum theories can be consistently defined by means 
of an adiabatically switched interaction, and to do it one 
has to allow time-dependent metrics to be used. Practical 
calculation of the scattering matrix does not necessarily 
entail solving (fT5"|k it is enough to calculate the static 
metric for the full Hamiltonian and then proceed with 
(|29| . Some complications may arise from the necessity 
of calculating the M0ller operators separately instead of 
the S-matrix as a whole; consequently one has to make 
necessary modifications to the standard methods of S- 
matrix evaluation. More importantly, one must be sure 
that the metric one has inserted into (|29|) is the same 
as the metric which arises from unity via the adiabatic 
switching. This presents an important restriction since 
it removes the ambiguity which was often encountered 
in the 'PT-symmetric theories. (The ambiguity can be 
removed by demanding simultaneous quasi-Hermiticity 
of an irreducible set of observables, but in field theory 
one has usually only the momentum and spin operators 
which commute with the Hamiltonian and so arc useless 
for the purpose.) 



We have argued that the necessary conditions for valid- 
ity of the described approach are met when the Hamilto- 
nian is diagonalisable and has finite gaps in the spectrum. 
We do not know to what extent the framework is appli- 
cable when continuous spectrum is present. Hopefully 
there is a consistent generalisation to reasonable class of 
models with continuous spectrum since there are versions 
of the adiabathic theorem that do not use the gap con- 
dition, see e.g. 
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